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1 Introduction 

In biomagnetic signal source activities, some are 
correlated, some are uncorrelated, and others are 
indifferent. If MEG (Magnetoencephalography) 
data have high signal to noise ratio, then the 
dominant components of the data are generated 
from such correlated, uncorrelated, and indifferent 
signal source activities. In this case, clustering of 
the signal sources in correlation and uncorrelation 
senses may be useful. Let the brain be discretized to 
small lattice points having current dipoles and a 
certain forward model for calculating magnetic 
field over the brain from a current dipole be 
assumed, the brain can be described mathematically 
by a linear system equation such that the clustering 
can be treated as simply grouping of the discretized 
elements. In this paper, we propose a concept of a 
biomagnetic inverse analysis with care for 
uncorrelated signal decomposition and for signal 
source localization simultaneously. 

The well-known related studies are mostly 
categorized in the following three: 

1. Uncorrelated / independent signal decom¬ 
position with no care for signal source 
localization; PCA (Principal Component 
Analysis), ICA (Independent Component 
Analysis), etc. 

2. Decomposition with no care for uncorrelation; 
beamformers (minimum variance beamformer, 
signal space projection, virtual electrode), etc. 

3. Signal source localization with no care for 
uncorrelation; equivalent current dipole 
method, minimum norm solution, SAM 
(Synthetic Aperture Magnetometory), MUSIC 
(Multiple Signal Classification), etc. 

2 Methods 

2.1 Model 

The brain is discretized spatially and the forward 
model employs the Sarvas equation [1], Therefore, 
the elements in each discretized point are a pair of 
linearly independent dipoles. The signal activity 
(magnitude of the vector) at the element p 


ip = 1 ,—,P) and at discrete time n (n = is 

denoted by s p [n] and its vector representation is s[«] 
= [ Sl [n],-,Sp[n]] J . 

The subset Q = {ca\n\, -, c, c Ui [n] } e s[w], 
(i = 1, •••,/, j = 1 ,■■•, J t ) consists of /th cluster. The 
elements in the same cluster are correlated with 
respect to each other, 

Vc, [«]e C,, c[«] = [cj[n], •,c / [«]] T 

=> ^c[«]c[«] T ^ ~ diag., 

where <•> is time average. The leadfield generated 
by the element p is denoted by \ p (time invariant) 
and its matrix representation is L = [li,—,lp]. The 
MEG time series of mth sensor (m = 1,—, M) is 
denoted by x m [n\ and its vector representation is 
x[n\ = [xi [«],—, x M [n\\'. As the result, the linear 
system equation is expressed, 

x[«] = Ls[«], (2) 

which is an usually under-determined equation such 
that the minimum norm solution exists. Since the 
x m [n] is highpass filtered, its time average <x m [n}> 
is almost zero. 

2.2 Beamforming 

Using the minimum norm solution of Eq. 2, the 
estimate s[«] of s[«] is calculated. For the stability 
of the solution, avoiding estimation for subspace 
corresponding to relatively small singular values of 
L is employed. However, the number of the 
avoidance should be greater than the effective rank 
Q x of <x[«]x[n] T >. 

2.3 Clustering 

Clustering itself has many methods. The aim of this 
paper is to propose the concept of the biomagnetic 
inverse analysis based on decorrelation. Therefore, 
the strategy of the clustering is only outlined as 
follows. First, the correlation coefficient matrix 
R = [/?J = [^[«].?>[«] 1 /1.?,[«] || |] of s\n | is 

calculated. Second, two thresholding values t\ (~1, 
for correlated) and t 2 (~0, for uncorrelated) are 
determined. Comparing the absolute values |/?*/|, 
find clusters having inter-correlation less than t 2 . 



Finally, eli mi nate elements in the clusters, which 
have intra-correlation less than t\. As the results, the 
elements of s[«] are divided into the clusters 
having high intra-correlation and low inter¬ 
correlation with respect to each other and the 
remaining indifferent cluster. If necessary, the 
indifferent cluster can be further divided under the 
same strategy only for the elements in the cluster. 
Furthermore, of course, the clusters consists of the 
discretized signal sources such that the localizations 
are also achieved simultaneously. 

3 Results 

3.1 Experiment 

Binaural AEF (Auditory Evoked Field) related to 
sound localization brought by ITD (Interaural Time 
Difference) was used for the evaluation as a sample. 
The stimuli for both ears consist of A (silent pair 
for both ears), B (correlated pseudo noise pairs 
having time-shift=0), C (same as B except time- 
shifteO), and D (uncorrelated pseudo noise pair). 
For example, in case of the stimulus B to C, the 
subject felt the movement of the sound source. The 
AEF was measured with a whole-head 64ch 
biomagnetometer system (CTF Systems Inc.). A 
typical stack plot and corresponding iso-contour 
map are shown in Fig. 1. 
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Figure 1: MEG related to the binaural AEF with 
ITD. a) Stack plot, b) Iso-contour map at the 
latency indicated by the arrow in a). 


3.2 Clusters of the signal source activities 

We applied the biomagnetic inverse analysis based 
on decorrelation to the AEF data. PCA analysis of 
x[n\ indicated that 4 dimensional subspace could 
be enough for covering over 90% power of the 


signals ( Q x ~ 4). Since the clustering was based on 
decorrelation, the number of clusters was also set to 
4. The stack plots and the 3D arrow maps 
corresponding to the clusters are shown in Fig. 2. 
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Figure 2: The biomagnetic inverse analysis based 
on decorrelation. Stack plots and 3D arrow maps of 
the clusters 1 to 4 are shown in a) to d) respectively. 


4 Discussion 

A typical AEF activity characterized by activities in 
both temporal lobes can be seen in Fig. 1. As 
shown in Fig. 2a, the clusters are distributed around 
temporal lobes. Therefore, the cluster 1 seems to be 
the dominant activity. However, the localization 
seems to suffer from typical minimum norm 
solution problems. Namely, the localized signal 
sources are apt to be distributed near the model 
surface and the distribution is not so focused. 
Although minimum variance beamformer was used 
instead of the minimum norm solution, the results 
were very sensitive to regularization methods. 
Signal decomposition methods like PCA and ICA 
usually have difficulty in interpretation of each 
component. On the other hand, beamformers cannot 
cluster the discretized signal sources and have also 
difficulty in the interpretation. Since this method 
derives not only the signal source activities 
(components) but also the locations of the clusters 
in correlation and uncorrelation senses, a priori 
physiological and anatomical knowledge might be 
directly helpful for the interpretation. 

This method will be beneficial to analysis for multi¬ 
modal activities in the brain. 
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